Maximal Violation of Bell Inequalities using Continuous Variables Measurements 
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We propose a whole family of physical states that yield a violation of the Bell CHSH inequality 
arbitrarily close to its maximum value, when using quadrature phase homodyne detection. This 
result is based on a new binning process called root binning, that is used to transform the continuous 
variables measurements into binary results needed for the tests of quantum mechanics versus local 
rcalistic theories. A physical process in order to produce such states is also suggested. The use of 
high-efficiency spacelike separated homodyne detections with these states and this binning process 
would result in a conclusive loophole-free test of quantum mechanics. 

PACS numbers: 03.65.Ud, 03.67.-a, 42.50.Dv 



I. INTRODUCTION 

Non-separability, or entanglement, has emerged as one 
of the most striking feature of quantum mechanics. In 
1935, it led Einstein, Podolsky and Rosen to suggest J|] 
that quantum mechanics is incomplete, on the premise 
that any physical theory of nature must be both "lo- 
cal" and "realistic". To quantify the debaté between 
quantum mechanics and local realistic (classical) theo- 
ries, Bell introduced a set of inequalities that must be 
obeyed by any local realistic theory whereas they are 
violated by quantum mechanics pi H, m. These results 
shifted the debaté from the realm of philosophy to expe- 
rimental physics. The experiments done at the beginning 
of the 1980s by Aspect and coworkers ||, |], |7| convin- 
cingly supported the predictions of quantum mechanics, 
but admittedly left open two so-called "loopholes" , that 
have to be addressed for the evidence to be fully conclu- 
sive. 

The first of these loopholes, called "locality"loopholc, 
arises when the separation between the measured states 
is not large enough to completely discard the exchange 
of subluminal signals during the measurements. The se- 
cond loophole, called "detection-efficiency"loophole, oc- 
curs when the particle detectors are inefficient enough so 
that the detected events may be unrepresentative of the 
whole ensemble. In 1998, Zeilinger et al. Q achieved 
communication-free condition by using a type-II para- 
metric down-conversion source and fast random switch- 
ing of the analyzers, that were separated by about 400m. 
This closed the "locality" loophole, but their detection ef- 
ficiency was not sufficient to close the second loophole. In 
2001, Rowe et al. Q measured quantum correlations bet- 
ween two entangled beryllium ions with up to 80% detec- 
tion efficiency, closing the "detection-efficiency" loophole, 
but unfortunately the ions were too close (about 3/íto) to 
avoid the "locality" loophole. Hcncc a present challcnge 
is to design and perform an experiment that closes both 
loopholes to lead to a full logically consistent test of any 
local realistic theory. 

Quantum òptics suggests good candidates, as pho- 



tons can be transported to sufficient long distances to 
avoid the "locality'Toophole. To close the "detection- 
efficiency" loophole, an alternative to photon-counting 
schemes consists in quadrature-phase homodyne mea- 
surements, that use strong local oscillators detected by 
highly efficient photodiodes. Up to date, a few theo- 
retical proposals that use quadrature-phase homodyne 
detections have been made |ï(], [ÏÏ], [Ï2| [Ï3| but for 
these set-ups the Bell inequality violation is a few per- 
cents only, that lies far away of the maximal violation 
attainable : 2y2 (compared to a classical maximum 
of 2) and (1 + \/2)/2 (compared to a classical maxi- 
mum of 1) respectively for the Clauser-Horne-Shimony- 
Holt (CHSH) (I and Clauser-Horne (CH)g inequalities. 
Gilchrist et al. |ï^, |ïï| use a circle or pair coherent state 
produced by non-degenerate parametric oscillation with 
the pump adiabatically removed. This state leads to a 
theoretical violation of about 1.015 (> 1) of the CH Bell 
inequality. Munro |Ï3] | considers correlated photon num- 
ber states of the form 
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n=0 
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where N is truncated at N = 10. He then performs a 
numerical optimization on each c n coefficient to maxi- 
mize the violation of the CHSH Bell inequality when an 
homodyne measurement is performed. For this specific 
state, the CHSH inequality is violated by 2.076 (> 2) 
and the CH inequality by 1.019 (> 1). From a different 
phase space approach, Auberson et al. fÏ4|| derive phase 
space Bell inequalities and propose a state that yields a 
maximal violation of up to 2\/2 (> 2). This state can be 
expressed in the position space by 
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where f(q) is a regularized form of 
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(2) 
with 



/_ dqf(q) 2 = 1. The main problem with this wave 
function lies in its singularities and phase switches. 
Therefore, it requires nontrivial regularization proce- 



II. ROOT BINNING 
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To begin our study, we consider a state of the form of 
a superposition of two two-particles wave functions with 
a relative phasc. 



i*> = ^(i//> 
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with / real, even and normalized to unity while g is as- 
sumed real, odd and normalized to unity. This kind of 
state looks similar to the one used by Auberson et al. 
]Ï4| , but the / and g functions will be quite different 
as well as the binning of the continuous variables being 
measured. 

The quadrature-phase homodyne measurement out- 
puts a continuous variable, and yet the majority of tests 
of local realism versus quantum mechanics require a bi- 
nary result. Hence for a given quadrature measurement 
qi (i = 1, 2) at either location A or B, we need to classify 
the result as either "+"or In ref. 0, [D| [u), the 
"positive-negative" binning is used, that is the result is 
classified as "+"if qi > and "-"if q t < 0. 

The choice of binning is quite arbitrary. For the state 
(||) we can consider another type of binning, we call root 
binning, that depends on the roots of the functions / and 
g (that are known in advance to the experimenters). We 
assign "+"when the result qi lies in an interval where 
f(q) and g{q) have the same sign, and "-"if qi is in an 
interval where / and g have opposite signs. We define 
D + as the union of the intervals in which f(q) and g(q) 
have the same sign and D~ as the union of the intervals 
in which /(ç) and g{q) have opposite signs. 

D+ = {V q e R | f(q) g(q) > 0} (6a) 
D' = {V q e R | f(q) g(q) < 0} (6b) 

Let us first consider the case when quadrature mea- 
surements in position space have been performed on both 
sides. So the binary probabilities we need for the CHSH 
type of Bell incquality will bc 
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FIG. 1: Schematic of a generalized Bell experiment. The 
source generates correlated states that are directed to the 
A and B devices used to perform the measurements, with 
adjustable parameters a and b. Each measurement provides 
a binary result "+"or "-" individually. 



with 



P{qi,q 2 ) = \{qi\{q 2 m\ 2 = 

\[f{qi?f{q2?+g{qi) 2 g{q2) 2 
+2 cos6 f(q 1 )g(q 1 )f(q 2 )g(q 2 )] 



(8) 
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Hereafter, we calculate the correlation function for (0): 



E„ 



P++ + P— - P+- - P-H 



(9) 



and as we have chosen / even and g odd, we get the 
remarkably simple expression : 



E quq2 =V 2 cose 



(10) 



where 



V 



f(q)g(q)dq- 

D + 

\.f(q)9(q)\dq 



f(q)g(q)dq 



D- 



(11) 



A similar binning will be applied for the momentum 
part. Since we suppose that f(q) is a real and even 
function whilc g(q) is real and odd, f(q) has a real 
even Fourier transform f(p) while g(q) has an imaginary 
Fourier transform ih(p), where h(p) is a real and odd 
function. Using these properties, and taking care of the 
supplementary i factor, the same reasoning applies for / 
and h as for /, g. Denoting as D' + and D'~ the intervals 
associated with / and h, we obtain 



-W 2 cosí 



(12) 



where 



W = I f(p)h(p)dp- / f(p)h(p)dp 

I D'+ JD>- 

\f(p)h(p)\dp (13) 



and equivalently, 



E quP2 = - V W sin( 



(14) 



E puq2 =-VW sm9 (15) 

Hence by combining Eqs.(^), ([Ï2]), ( fÏ4"|) and ( |Ï5| ) we can 
write the CHSH inequality(|) 



| cos(6)(V 2 + W 2 ) - 2 sm(0)VW\ < 2. 



(16) 



The maximum of S with respect to 9 is obtained for 
tan0 m = -2VW/{V 2 + W 2 ), and we have m -> -tt/4 
as V, IF — > 1. For this optimized # m we get the Bell 
inequality 



S = WW i + V 4 + 6V 2 W 2 \ < 2 



(17) 



Using this really simple expression ([H]), the debaté bet- 
ween quantum mechanics and local realistic thcories boils 
down to find functions / and g whose integrals V, W 
violate (^7|). An interesting feature appears when the 
distributions are eigenstates of the Fourier transform, so 
that V = W and equation (fT?]) becomes 



S = 2V2 V 2 < 2 



(18) 



So if such functions have the right overlap needed to ob- 
tain V = 1, one will get the maximal violation of the 
above inequality, which is obtained for S — 2v / 2· 

Whcn compared to the positive-negative binning, root 
binning has the advantage of having two parameters V 
and W to play with while the positive-negative binning 
has only one |Q. Moreover, as we will show now, the 
above Bell inequality is violated by simple wave func- 
tions, that no longer have the singularities that appeared 
inref.Q. 



III. FOUR PAWS SCHRODINGER CATS 

In order to propose an explicit expression of a state 
that violates the Bell inequality (|Ï7|), let us first consider 
the case of a Schròdinger cat state. The Schròdinger cats, 
that are a superposition of two coherent states of ampli- 
tudes a and — a, involve intrinsic quantum features such 
as negative Wigner functions, which make them interest- 
ing candidates for our case (||). We must choose an even 
cat for f(q) and an odd cat for g(q) : 

/( g )(xe-(«+") 2 / 2 +e-('- a ' 2 / 2 (19a) 
g(q) cx _ e -(«+ a ) 2 /2 + e -(i-) 2 n (i 9b ) 

unfortunately for this simple state we get V = 1 and W ~ 
0.64 for a -> 00, so that S ~ 1.90 < 2. Therefore this 
state cannot be used for violating Bell inequality (note 
that Gilchrist et al. |ï(], |ïï) also consider Schròdinger 
cats, but without getting a violation of Bell inequalities). 

Instead of Schròdinger cats which can be alive or dead, 
we consider acrobat cats which have 4 paws. Let us for 
instance consider: 

f(q) OC _ e -(?+3-) 2 /2 + e -( ? +a) 2 /2 + e -( q -af/2 

_ e -( q -3af/2 (2Qa) 



g{q) oc -e 



-(g+3a) 2 /2 e ~(q+a) 2 /2 + e -(q-af/2 
, p -(9-3a) 2 /2 



(20b) 



/(<?) and g(q) are depicted in Fig.^j together with their 
Fourier Transform. Note that for this choice of functions, 
each peak is distant from its neighbours of a = 2a, this 
disposition yields an optimal overlap of / and g and thus 
a high value of S. The best violation appears when the 
peaks move off as a — > 00. In that case, V = 1, W = ^ 
and thus we get the significant violation of S — 2 ~ 0.417 
(in facts, the condition a — > 00 appeared to be not so 
strict numerically, as an amplitude a = 5 is enough to 
obtain S ~ 2.417). Such a violation represents a large 
improvement compared to Munro's best result of 0.076 
|Ï3| , for a state with no singularity and at least as easy 
to produce as Munro's c„ optimized state (|ï|). However, 
we are still away from the maximal value 2^/2 of the 
CHSH Bell inequality ]Ï5| |. In the following section, wc 
will propose a set of states to get closer to the maximal 
violation. 



FIG. 3: Infinite-paws states, represented in the position and 
momentum phase space. Thick segments denote Dirac Delta 
functions. 



goo overlap perfectly, these distributions would yield the 
maximal violation of the CHSH inequality, as S — 2\[2 
from eq.(p"8|). Of course, they are unphysical (non- 
normalizable) states, but the wave functions /jv and g^ 
in (|2Ï|) can be considered as regularized forms of and 
goa, widening the Dirac Delta functions to Gaussians and 
taking a finite number of paws. Thus one can understand 
that S -> 2^2 as N -> oo. 

Another regularization of the wave functions ( p2| ) con- 
sists in a widening of the Dirac functions to Gaussians 
of width s associated with a Gaussian envelope of width 
1/s: 

foo- a ,s(q) oc Gy a (q) [foo-a * G s (q)} (23a) 
(q) oc G\ /s (q)[goc-, a * G s (q)] (23b) 

2 

where G s (q) — exp(— jt), * denotes the convolution 
product, and s is a squeezing parameter. When s <C 1, 
one indeed has: 

foo:a,s(p) (XG S * [/^/«Gi/J (p) « /oo;7r/a, S (p) (24a) 
s (p) OC G s * [flf o;ir/oG ! l/«](p) ~ 9oo;tt /a,s ÍP) (24b) 

The violation exhibits a symmetry, as 5(a) = S'(— ), with 
a maximum reached for a = y 7 ^? when / and g are ap- 
proximately eigenstates of the Fourier transform. 

Thanks to the Gaussian envelope, the above functions 
can be truncated to a finite total number of paws N with- 
out affecting numerically the / and g functions, provided: 

2J2\ lnel , x 

JV > -X— ■ [ 25 

as 

where e is an arbitrary small tolerance parameter. For 
e = 0.01, a = ^/ïr and s — 0.3, the condition yields 
N > 12. Given these parameters and iV = 12, we get 
S ~ 2\/2 with a relative error of 0.01%. Regardless to 
this condition, one may also arbitrarily choose to limit 
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The preparation procedure begins with the quantum 
system in the vacuum state |0) and an ancilla qubit in 
the ground state |0) a (if some squeezing parameter is 
needed, one may take a squeezed vacuum as quantum 
system and follow the procedure we describe here. For 
simplicity reasons, we set s — 1). Let us first apply: 



H e 



H 



(26) 



with p is the momentum operator applied to the contin- 
uous variable state, H the Hadamard gate and o~ z the 
Pauli matrix applied to the qubit : 




1 

-1 



(27) 



We then have a probability 1/2 of measuring the qubit 
either in the excited or in the ground state. If it is found 
in the \l) a state, the continuous variable is left in the 
state |Ti) oc — | — a) + \a), otherwise the procedure is 
stopped and we try again. The qubit is then bit-flipped 
to |0) a and we go on by applying the sequence 



H e 



-i2apa z 



H 



(28) 



Measuring the qubit in the |0) a results in the continuous 
variable left in the state IT2) oc — | — 3a) + | — a) — \a) + 
\3a). To increase the number of paws, we iterate the 
following procedure given |T ra _i) and the qubit in |0) a : 



Apply the operators 



H e 



-i2™~ 



H 



(29) 



FIG. 4: N — 12 cat states described by eqs(p3|) presented 
in position space (left) and in momentum space (right) for 
parameters a = 1.8 and s — 0.3. Left axes are in arbitrary 
units and normalized to unity. 



V. CONDITIONAL PREPARATION OF 
ENTANGLED N PAWS SCHRÒDINGER CATS 

Preparing entangled TV paws Schròdinger cats is a 
challenging task, we begin by focusing our attention on 
how one single N paws Schròdinger cat as defined by 
(211) could be generated. This state has strong similari- 
ties with the encoded states introduced by Gottesman 
and coworkers |ïl| to perform quantum error-correction 
codes. Recently, Travaglione and Milburn jÏ7| presented 
a proposal to generate non-deterministically such encoded 
states. Following this study, we will first show how to 
generate the state \g) by applying a specific sequence of 
operations similar to |l7| and then derive a setup to pro- 
duce the whole state (Eiï. 



• Measure the qubit 

• If the qubit was in the state |0) a , we have created 
|T„) 

• If the qubit was in the state |l) a , discard and try 
again. 

Once the number of paws is considered satisfactory, we 
stop the previous iteration. The last point to generate 
\g) is to apply the following sequence 



H e 



H 



(30) 



If the qubit is found in the state | 0) a , we have cre- 
ated the state \g) as defined in eq.( 21b ). If the pro- 
tocol was stopped after n iterations (corresponding to 
|ï„)), the state \g) is created with probability l/2" +1 
and shows N = 2 n+1 paws. In the last section of ref . , 
Travaglione and Milburn briefly consider the question of 
the physical implemcntation of this iteration process us- 
ing a radio- frequency ion trap. We refer the reader to 
this article for further details. 

For the creation process of the state |ï') = 
"^Tf (I//) + ete \99))i we present on fig.|| a global view 
of our scheme. The T-labelled box corresponds to the 



6 



|ff>+e ie |gg> 




FIG. 5: Schematic of the setup used to generate the state 
l*> = ^5 (|//> + é e \gg)) with / and g defined by eqs.©. 
See text for the notations. 



generation of \g) we just saw. Our process is based 
on a CNOT gate that entangles two qubits {\ip) a — 
-^(\íl) a + e %e \00) a)) ■ Each qubit is then associated with 
a state \g) through the following operator 



A = 



(31) 



where q is the position operator of the \g) state and a z 
the Pauli matrix applied to the qubit. If the qubit is at 
zero, the \g) state will be left unchanged. If the qubit is 
at one, for s <C a the sign of each two peak is changed 
so that the \g) state is changed into the |/) state. After 
this operator A has been applied, the whole state is left 



\g g 0) + |/ / 1 1) 



(32) 



In order to disentangle the two qubits from the conti- 
nuous part, one has to make the qubits pass through 
an Hadamard plus a CNOT gate. Finally measuring one 
qubit in the zero state will project the continuous variable 
system onto the awaited state 



l*> = ^(l//> 



' i9 \99)) 



(33) 



Such a process generates the N-paws Schròdinger cats 
as defined in (H|) where each peak has the same height. 
The case (|2^) with a gaussian envelope is more compli- 
cated to produce because each peak must be separately 
weighted, but as seen in table |, the violation for a state 
( pT| ) with equal heights is already quite strong (S ~ 2.41 
foriV = 4). 



VI. CONCLUSION 

Considering quadrature phase homodyne detection, we 
have derived a new binning process called root binning to 
transform the continuous variables measured into binary 
results to be used in the test of quantum mechanics ver- 
sus local realistic theories. For this process, we propose a 
wholc family of physical states that yield a violation ar- 
bitrarily close to the maximal violation in quantum me- 
chanics and much stronger than previous works in the 
domain 0, |Ü|, 0, g|. 



We have also tested root binning on other interest- 
ing forms of Bell inequality that are the information- 
theoretic inequalities dcvclopcd by Braunstein and Caves 
|Ï8| and generalized by Cerf and Adami |ïs|j . Using 
quadrature measurements, root binning and our state 
defined by eqs.(^|p^), we unfortunately could find no 
violation for neither Braunstein's nor Cerf's form of 
information-theoretic Bell inequality. Our state in fact 
tends to the minimum limit for violation of these informa- 
tion inequalities when V — ► 1 and W — ► 1. As a matter 
of fact, the binning process discards a lot of information 
that lies within each interval of the binning. This in- 
formation loss may prevent any violation of information- 
theoretic inequalities. 

As a conclusion, let us point out that though the 
present idea sounds quite attractive, its practical implc- 
mentation is very far-fetched. Though we do propose a 
theoretical scheme to prepare the required states, it re- 
lies on various features (coupling Hamiltonian, CNOT 
gates) that are not presently available with the required 
degree of emciency. For going in more details into the 
implcmentation of the present scheme, inefficiencies and 
associated decoherence effects should be examined in dc- 
tail for each required step, i.e. preparation, propagation, 
and detection. Also, various possible implementations of 
the proposed scheme should be considered |)| [Ï7], ^ï) . 

Such a study is out of the scope of the present paper, 
that has mostly the goal to show that arbitrarily high 
violations of Bell inequalities are in principle possible, by 
using continuous-variables measurements and physically 
meaningful - though hardly feasible - quantum states. 

This work is supported by the European 
IST/FET/QIPC program, and by the French pro- 
grams "ACI Photonique" and "ASTRE". 



APPENDIX A: EXPRESSION IN THE FOCK 
BASIS 

An interesting point is to study the decomposition of 
states |/) and \g) on the Fock basis |n), which are ex- 
perimentally accessible. Starting from states (^3j) with 
a = t/tt, s = 0.4 and with N satisfying condition (|2£ 
we obtained after truncating at the 14 th order: 

|/) = V0.459|0) - V0.491|4) - V0.008|8) 



- V0.042|12) (Ala) 
\g) = V0.729|l) + V0.155|5) - V0.107|9) 

- V0.009|13) (Alb) 

These states allow to reach a violation S — 2 = 0.81. The 
state |/) only involves orders n = (mod 4), as we have 
n = l (mod 4) for \g). This directly comes from the fact 
that (p\n) = (—i) n (q\n) q=p and (—i) 4 = 1, so that states 
of the form J2 n =a (mod 4) c «l n ) are eigenvectors of the 
Fourier transform. From these considerations we have 
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obtained S — 2.68 for: These states are quite simple, and it is expected that a 

specific procedure to producc thcm might be designed. 

|/) = V0.585|0) - V0.415|4) (A2a) 
\g) = V0.848|l) + V0.152|5) (A2b) 

and S = 2.3 for: 

= VrT67|0)-VO33|4), |.g) = |l) (A3) 
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